Thermoelectric transport of mesoscopic conductors coupled to voltage and thermal 

probes 
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We investigate basic properties of the thermopower (Seebeck coefficient) of phase-coherent conduc- 
tors under the influence of dephasing and inelastic processes. Transport across the system is caused 
by a voltage bias or a thermal gradient applied between two terminals. Inelastic scattering is mod- 
eled with the aid of an additional probe acting as an ideal potentiometer and thermometer. We find 
that inelastic scattering reduces the conductor's thermopower and, more unexpectedly, generates a 
magnetic-field asymmetry in the Seebeck coefficient. The latter effect is shown to be a higher-order 
effect in the Sommerfeld expansion. We discuss our result using two illustrative examples. First, 
we consider a generic mesoscopic system described within random matrix theory and demonstrate 
that thermopower fluctuations disappear quickly as the number of probe modes increases. Second, 
the asymmetry is explicitly calculated in the quantum limit of a ballistic microjunction. We find 
that asymmetric scattering strongly enhances the effect and discuss its dependence on temperature 
and Fermi energy. 

PACS numbers: 73.23.-b, 73.50.Lw, 73.63. Kv 



Introduction. Recent advances in heat measurements 
have enabled to envisage promising thermoelectric ap- 
plications at the mesoscale."'^ In particular, there have 
been a number of exciting proposals ranging from ther- 
mal analogs of electronic rectifiers.^ and transistors^ to 
efiicient converters of heat into electricity'* which can be 
shown to reach optimized configurations.^ Of fundamen- 
tal importance is the experimental verification of the car- 
riers' charge sign using thermopower techniques onlyj^ 

Since these novel devices operate at the quantum 
regime of transport, it is of primary interest to investigate 
dephasing effects which may be detrimental to their per- 
formance. Additionally, energy flow at finite temperature 
is expected to be particularly sensitive to inelastic transi- 
tions inside the sample due, e.g., to coupling to phonons. 
This problem has been addressed only very recently j^ 

A convenient way to introduce dephasing and energy 
relaxation in a mesoscopic conductor is based on the volt- 
age probe model.^ A fictitious terminal is attached to 
the sample such that the net current fiowing through the 
probe is set to zero. Hence, a carrier absorbed by the 
probe is reemitted into the conductor with an unrelated 
phase. The clear advantage of this approach lies on its 
simplicity and its independence of the microscopic details 
of the phase-randomizing processes. Therefore, universal 
properties of generic conductors can be investigated this 
wayjiSiiii The model has been recently applied to find the 
temperature and chemical potential profiles of an array 
of quantum dots J^ 

Including a third terminal acting as a dephasing probe 
requires to extend the scattering approach to the mul- 
titerminal case. For thermoelectric transport, this was 
achieved by Butcher. ^^ Experimentally, there are already 
available data for paradigmatic mesoscopic systems such 
as two-terminal point contactsi^ and chaotic cavitiesj^ 
The multiterminal case, however, has been less explored. 



A very recent exception is Ref. Il6| where a voltage drop is 
generated transversally to a horizontal thermal gradient 
in a four-terminal setup with an asymmetric scatterer. 

Here we investigate incoherent scattering effects on the 
thermopower (Seebeck coefficient S), defined at linear 
response as the ratio of voltage AV to temperature A9 
differences applied between two terminals (1 and 2 in Fig. 
[1]) at vanishing current /, 
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in the presence of an attached probe (terminal 3 in Fig. 
[T]) acting as an ideal potentiometer and thermometer. 
Quite generally, we find that S decreases as compared to 
the case without probe and, strikingly, the presence of 
incoherent scattering due to the probe causes the devel- 
opment of a magnetic-field asymmetry in S. To illustrate 
our findings, we analyze (i) a mesoscopic conductor with 
generic properties (a chaotic cavity); and (ii) a ballistic 
wire in the quantum limit (a few propagating modes) 
with an asymmetric scatterer. In the latter case, we 
explicitly calculate the size of the magnetoasymmetry. 
This result is relevant in view of recent predictions that 
relate this asymmetry to efficiency limits of thermal de- 
vices in converting a temperature gradient into electrical 
workJIi^ 

Model. We consider a mesoscopic conductor coupled to 
three terminals. Thermoelectric transport is described by 
the electric current, /„, and the heat current, Qa, flowing 
from terminal a = 1,2,3. In the linear response regime. 
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FIG. 1: Sketches of (a) a chaotic cavity attached to voltage 
and thermal reservoirs and (b) a microjunction with an asym- 
metric scatterer coupled to a probe as a model of a quantum 
wire with incoherent scattering. 



the transport equations read: 

la = 2^ GapVfj + 2_^ LapOp , 



(2) 

(3) 



where eVa = Ep — fj,a are voltage shifts away from the 
common Fermi energy Ep, with fia the electrochemical 
potential of lead a, while 9a = 0—Ta measures deviations 
of the reservoir temperature T^ from the common (bath) 
temperature 0. The transport coefficients are expressed 
in terms of the transmission probability Tap'. 



Gap = ~gv / dE {Na Sap - Tap)!^ ., 
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Kaf3 = g fdE (E - Ei 



f{Na5ap-Tap)f'o, (7) 



with gv — 2e^/ft. the conductance quantum, Na the num- 
ber of propagating channels in lead a and /q the energy 
derivative of the Fermi distribution function evaluated 
at Vq. = 6'q. = 0. In Eq. ([6]) we have used the Onsager 
reciprocity relations between cross terms. 

Isothermal case. We first consider the case where the 
temperature probe is an externally fixed parameter. This 
would correspond, e.g., to a phonon subsystem main- 
tained at a different temperature.^ The probe thus works 
as an ideal voltmeter with voltage V3 determined from the 
condition /a = 0, 
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We substitute Eq. ^ into the expression for /i, from 
which we find the current I ~ Ii ~ ~l2 flowing through 
the system: 



-I — (]Vi2+Li29i2+Lij,9ii- 



Gi 



G31 + G32 



-[^31^13+^32^23] • 
(9) 



7 is a function of voltage Vij = Vi — Vj (i, j = {1, 2,3}) 
and temperature 9ij = 9i — 9j differences, as should be. 
In Eq. di]), -g = G12 + Gi3G32/(G3i + G32) is the well- 
known expression for the conductance in the presence of 
a voltage probe in the purely electric case for which all 
temperature shifts are set to zerOi^ 

In the isothermal configuration, the probe is main- 
tained at the same temperature as the common bath 
(^3 — 0). We are then free to specify the precise form 
of the temperature gradient 9i — 92 y^ 0. We choose the 
symmetric bias 9i = — ^2 = ^9/2, which is commonly 
employed in actual measurements J^ We then compute 
the thermopower from Eq. ([IJ for V\ — V2= AF: 
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Equation (ITU|) shows two contributions to the ther- 
mopower as compared to the case without phase ran- 
domization for which S — — ii2/Gi2. The first part cor- 
responds to direct inelastic scattering with the probe and 
is proportional to the term L13, as can be expected from 
an analogy with the purely electric case (cf., g above). 
The second part is more surprising — it is nonzero only 
if there is an asymmetry between the transmissions into 
the probe of those carriers injected from lead 1 and lead 
2 (^32 7^ ^31 )■ This term has no counterpart in the 
purely electric case. A similar asymmetry effect, but 
with phonon carriers, has been recently shown to be cru- 
cial in the development of rectification effects in dielectric 
junctions. ^^ 

Adiabatic case. More interesting is the case when the 
probe plays simultaneously the role of an ideal poten- 
tiometer and thermometer. Then, not only the charge 
current but also the energy flux is zero at the probe. 
These two conditions determine V3 and 6*3.'^" Recently, 
it has been suggested that electronic local temperatures 
can be consistently defined introducing thermal probcs3i 

Imposing 73 = and Qt, = 0, wc find: 
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which we substitute in the equation 7i = in order to 
obtain the thermopower S = AV/ A9 = 1^12/^12, 

S = [Gi3(L3i7ir33 — 7ir3iL33) + G33(Li37f3i — Lii7ir33) 
+^-^33(-^13-^31 ~ -^11-^33)] / [^-^33(Giii33 — G31L13) 
+ G33(6'77i3L3i -I- Gii74r33) — Gi3(G3i7ir33 -t- 6'L3iL33)] . 

(13) 

Chaotic cavity. We now consider a generic mesoscopic 
sample — a metallic quantum dot whose classical analog 
displays chaotic dynamics [see Fig. UKa)]. Its isotropic 



properties permits us to treat the cavity as an effectively 
zero-dimensional object characterized with a mean level 
spacing S. Transport occurs when the cavity is coupled 
to external reservoirs usually via quantum point con- 
tacts with a large number (iVi and N2) of propagating 
channels. The experimentally relevant case deals with 
clean samples.—^ As a consequence, transport is ballistic 
and its corresponding statistics is well described by ran- 
dom matrix theory (RMT)i^ In what follows, we assume 
eVa, ksOa < Et, with Et = {Ni+N2)5 the Thouless en- 
ergy. 

Using a Sommerfeld expansion, one finds approximate 
expressions for the response coefficients in Eqs. (|4])-([7]) 
in terms of the transmission functions T^p and their 
energy derivatives T'^a evaluated at the Fermi level: 
Gap — 9v{Na5ap — Tajs), Lafj = QoT^p and Kap — 
-geiNJai3-Taf3)/e,wheTegg = (2e//i)(7r2fc|6l/3). Note 
that the ratio 0gv/egg = 7r^fc^/3 depends on universal 
constants only, representing the fundamental quantum 
of thermal conductance of a perfectly transmitting mode 
with averaged reservoir temperature 0.^'* 

We investigate the statistical properties of the ther- 
mopower, in particular its mean value and variance. 
In the limit Na ^ 1 we can calculate correlations 
of Tap and T^p within RMT. We consider the cases 
/3 — 1 (orthogonal ensemble) and /3 — 2 (unitary en- 
semble) corresponding, respectively, to the presence and 
the absence of time reversal symmetry. We define the 
functions A^^ = N^N^/Nf - S^^^N^/Nt and Bp„ = 
\/2{NpNa/Nf - SpaNp/N^), where ^, v, p, a are lead in- 
dices and Nt is the total number of modes. Then, using 
Ref. u3 we derive the useful relations. 
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From Eq. (J15l) we immediately see that the mean 
vanishes, (S) — 0. This result was found in Ref. i26 for 
the case without dephasing. The reason is clear: the 
L coefficients are functions of the energy derivative of 
the scattering matrix, but this is zero on average (the 
derivative can be positive or negative for a specific sample 
of the ensemble but fiuctuates on average around zero). 
The fluctuations are described by varS* — (S^) — {S)^, 
which are generally nonzero. For definiteness, we set 
Ni = N2 — N. Then, the thermopower fluctuations 
for both the isothermal and adiabatic cases are governed 
by the same expression: 



var S" 



87T^k%9'^ 1 + Sfii 
9e2(52 {2N + iVa)" 
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Our result generalizes the thermopower fluctuations to 
the case of finite dephasing. Note that the fluctuations 
are not universal and vanish quickly as the mode num- 
ber increases. Strikingly, the same functional dependence 



{N~'^) appears in the magnetoasymmetry of the weakly 
nonlinear conductance^! since it also depends on the en- 
ergy derivative of the transmission. The fluctuations 
are twice larger in the absence of magnetic fleld. Im- 
portantly, incoherent scattering effects reduce drastically 
the fluctuations of S since these have a quantum origin 
and the probe introduces decoherence in the system. To 
leading order, one has (varS* — var 5^^^")/ var S'^^^'^ — 
1-2N3/N + 0{N3)\ 

In the adiabatic case, we have neglected in S terms 
that do not contribute to the variance (e.g., terms that 
involve three L's are of higher order). Then, Eq. ([T^ can 
be approximated as. 
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Since for symmetric couplings (A^i — N2 = N) the pref- 
actors of Lap in Eq. P^ are, on the ensemble average, 
equal to those of the isothermal case [Eq. (ITUl) ]. we obtain 
the same expression for var S. Thus, only asymmetric 
couplings (A^i 7^ N2) could distinguish between the two 
cases. 

Magnetic-field asymmetry. In a two-terminal con- 
ductor, the linear conductance is an even function of 
the applied magnetic fleld B. This fundamental sym- 
metry remains unchanged even after elimination of the 
probe coupled to a mesoscopic conductor since g{B) = 
g{—B). This can be seen by recasting g as g = Gn — 
G13G31/G33, which is manifestly symmetric under B re- 
versal. In contrast, while the two-terminal thermopower 
S = — Lii/Gii is manifestly an even function of B, this 
statement does not hold in the presence of a probe in the 
adiabatic casein. From Eq. (fT5|) we find the thermopower 
magnetoasymmetry $ = S{B) — S{—B): 

$ = [G33(Li3iir3i — L3iA'i3) + i33(G3ii4ri3 — Gi3A'3i) 
+^33(Gi3L3i — G3iii3)] / [6'L33(GiiL33 — G31L13) 
+ G33(6'Li3L3i + G11A33) — Gi3(G3iAr33 -I- 6'L3iL33)] , 

(20) 

which is generally nonzero. Then, incoherent asymmetry 
changes the symmetry of the Seebeck coefficient. 

We recall that the leading order in a Sommerfeld ex- 
pansion reads Kap — —{■7T^kj^9/3e'^)Gap- Substitut- 
ing in Eq. (PO)) . one finds $ = 0. As a consequence, 
-B-asymmetries in the thermopower are a higher- order 
effect. This implies that the size of $ will decrease 
quickly as temperature decreases (typically, as 9^).'^-^ The 
leading-order Sommerfeld approximation neglects terms 
of order {ksd/EpY^ This factor is rather small in 
macroscopic metals and it is thus not necessary to con- 
sider higher orders. But in low dimensional systems the 
B-asymmetry of S should be visible since Ep and 9 can 
be of the same order. We note in passing that the con- 
ductance g is, in contrast, always B-symmetric indepen- 
dently of the Sommerfeld approximation. 

Numerical simulations. Clearly, for a chaotic cavity 
the thermopower magnetoasymmetry vanishes on aver- 
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FIG. 2: (Color online) Numerical results for the thermopower 
(a,b) and thermopower _B-asymmetry (c,d) of a microjunc- 
tion with an asymmetric scatterer [Fig. IHb)] as a function of 
Fermi energy (a,c) and temperature (b,d). Full (open) sym- 
bols correspond to the presence (absence) of the voltage and 
thermal probe. Solid (dashed) lines in (b) are obtained us- 
ing the Sommerfeld expansion to lowest order with (without) 
the probe. Note the strong deviation at moderately low tem- 
peratures. Solid line in (d) shows the $ oc 0'^ dependence of 
the B-asymmetry at low temperature. Parameters: magnetic 
length hc/eB — L^, temperature ksO = f? /mL^ (a,c) and 
Fermi energy Ef = IQh^ /mL^ (b,d). 



age. Moreover, we expect the fluctuations to be exceed- 
ingly small since $ depends, to leading order in 0, on 
the second derivative of the transmission. Therefore, it is 
more natural to estimate the size of $ in a different meso- 
scopic system — a microjunction as in Fig. [IJb). Ther- 
mopowers yielding 0.6 /xV/K have been recently detected 
in a similar setupJ^ The asymmetric scatterer deflects 
the electronic trajectories differently depending on the B 
direction. Therefore, we expect an asymmetry in 5 when 
B is inverted. 



Our prediction is confirmed by the numerical calcu- 
lations of Fig. [21 We compute TapiE) using a grid 
discretization of the Schrodinger equation and, subse- 
quently, the response coefhcients from Eqs. (2])- ([7]) with 
Gauss-Legendre quadratures.^^ Upper panels show that 
adding the voltage and thermal probe reduces the ther- 
mopower absolute value as a function of both Ep and 9, 
in agreement with our previous results. In Fig. [DJb) we 
show that S deviates from the lowest-order Sommerfeld 
expansion as temperature increases both with and with- 
out the probe. We emphasize that the thermopower is 
_B-asymmetric only if the probe is present, as shown in 
Figs. [2lc,d). In Fig. [2]i) we plot the asymmetry $ as a 
function of temperature and show that $ increases as 9^ 
at low temperature. The maximum value of $ depends 
on the details of the scatterer and the system's parame- 
ters. We find that $ can reach values as high as 2% for 
Fermi energies close to the activation of the second mode 
Ef ~ 19.7ff/mL^. 

Conclusions. To summarize, we have investigated inco- 
herent scattering effects on the thermoelectric transport 
of a mesoscopic conductor using a fictitious probe acting 
as an ideal potentiometer and thermometer. Our main 
findings are: (i) a general expression for the quantum 
fiuctuations of the Seebeck coefficient upon elimination of 
the probe valid for both isothermal and adiabatic probes; 
(ii) a magnetic-field asymmetry of the thermopower, re- 
quiring both inelastic and dephasing processes, which be- 
comes apparent only when higher-order terms are consid- 
ered in a Sommerfeld expansion. Odd-in-J5 thermopow- 
ers in Andreev interferometers have been experimentally 
observed^ and theoretically studied)^ We hope that our 
results will motivate the experimental detection of asym- 
metric thcrmopowers in normal systems. 
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SUPPLEMENTARY MATERIAL 



Details on Numerical Calculations 



This section discusses in detail our numerical simulations for a few-mode mesoscopic conductor subject to incoherent 
scattering. We determine the set of multiterminal transmissions Tap{E) using a grid discretization of the Schrodinger 
equation and solving the resulting linear equations, varying the boundary conditions for incidence from each of the 
leads. The leads are modelled as narrow waveguides with impenetrable hard walls, forming a T-junction geometry as 
sketched in Fig. [TJb) of the main paper. We also include an asymmetric triangular scatterer located at the center of 
the junction. 

Magnetic field effects are modelled as follows. To avoid extended vector-potential fields A{f) inside the leads, we 
consider the localized fields of Fig. |31 which might correspond to two coaxial coils of radii i?i and i?2, perpendicular 
to the junction and carrying opposite currents. In this geometry B{f) = B{r)uz and A{r) = A(r)u^, where Uz and 
W0 are the Cartesian and the azimuthal unit vectors, respectively, while r is the distance from the junction center 
(Fig. 131). The field strengths read 



r <Ri 
Ri <r < R2 

r> R2 



B{t) 
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(21) 



Notice that, as required, not only B{r), but also A{r), vanishes in the leads when r > R2. Therefore, our results arise 
purely from magnetic effects that occur within the junction only. 




lead 1 



lead 2 



FIG. 3: (Color online) Sketch of the magnetic field modelling, with the two regions of radii Ri and i?2. 



The central scatterer is modelled as a smooth two-dimensional barrier of height Vq with a triangular shape in the 
plane. More specifically, the scatterer's potential barrier in polar coordinates is given by 



V{r,e) = 



Vn 



(22) 



1 + e('-^i'W)/'^T ' 

where Rt{0) is a piecewise function characterizing a triangle of side Lt and oriented along a specific angle 0t, and 
ctt a small diffusivity. 
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FIG. 4: (Color online) Density distribution, |^(r)| , with ^{f) the solution of the Schrodinger equation for incidence from 
terminal 1 (left). The density is given in units of L~^ , where L is the lateral extension of the leads. The position of the 
leads and the scatterer is indicated by solid lines. The incident energy is E = 30h^ /mL^. Parameters of the magnetic field: 
hc/eB — L^ , Ri = L, R2 = 1.5L. Parameters of the repulsive scatterer: Vb — 200ft^/mL^, Lt ~ 0.51/, ar ~ 0.05L, 6t = 15°. 



The numerical calculation of the kinetic terms including vector potential contributions is made using the Peierls 
substitution for lattice calculations. This is very convenient since it fulfills gauge invariance by construction. The 
practical implementation is as follows. We denote with tij the finite-difference matrix associated with the kinetic 



energy on the grid. Thus, 
neighboring points j: 



the operator acting on function 4* at grid point i is given by the following sum over 
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The Peierls substitution for the vector potential terms amounts to including the following field-dependent phases in 
the sum 
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2m 
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(24) 



As an example, Fig. |4] displays the probability density in the presence of scatterer and magnetic field in the case of 
incidence from the left (terminal 1). 

Once the wave function for each incidence condition is known, it is an straightforward task to determine the 
quantum transmissions Tq^ to be used in the evaluation of the linear transport coeSicients, Eqs. (4)-(7) of the paper. 
These integrals are performed using a controlled numerical method, with the required precision (10^^) achieved by 
successive partitioning of the integration domain in an increasing number of zones and computing each of them with a 
Gauss-Legendre method. In practice, the transmissions Tap{E) for arbitrary energy E are obtained by interpolating 
or extrapolating on a table of closely spaced precomputed values. Special care is taken for E values close to the 
transverse mode energies to keep the discontinuities at these energies of the transmissions. 



LOW TEMPERATURE DEPENDENCE 



In this section, we show that at low temperature the thermopower asymmetry has a 9^ dependence. We present our 
analysis, both numerical and analytical, of this limit. We recall that the thermopower _B-asymmetry $ = S{B)—S{—B) 



$ = 



reads, 

GssiLlsKsi — L31K13) + L33(G3iifi3 — G13K31) + K33{Gi3L3i — G31L13) 

^-^33(Gii£33 — G31L13) + G33(0Li3L3i + GwK-^'i) — GniGzxK-^'^ + 6'L3iL33) 
where the transport coefRcients are given by 

Gq^ = gv{NaSal3 — Tali) , 

Lap = geT'aji , 

Kali = -geiNaSap - Tali) / e , 



(25) 



(26) 
(27) 
(28) 



to leading order in the Sommerfeld expansion. Higher orders are 0{d)^. It follows from Eqs. ([26]) and psp that Kais — 
— {TT^kg/3e^)6Gaii- Substituting this result in Eq. (I^H)) one obtains $ = 0. Therefore, nonzero magnetoasymmetries 
arise only if we carry out a Sommerfeld expansion for Kap to higher order. We find. 



K, 



Q/3 



-^eG 



a/3 



77r%_^3j..^^ 0(^)5 



15/i 



'-aP 



(29) 



Since Lap is 0{d), the numerator in the right-hand side of Eq. (|25p is 0{9)'^, whereas the denominator is 0{9). Then, 
the leading order of $ is 0{9)^. 

Our numerical results are shown in Fig. [SJ which displays a fit of $ versus 9 (dots) to a cubic-law behavior (line). 
The behavior predicted above is in agreement with the numerical data at low temperature. Obviously, higher-order 
contributions start to dominate as temperature increases. For the asymmetric microjunction, we estimate that the 
thermopower asymmetry deviates from the 9^ law at fcs^ > h^ /mLp' . 
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FIG. 5: (Color online) Low temperature fit of the thermopower asymmetry, <& oc S'^. Logarithmic scale is used in the main 
plot to highlight the low temperature region. For comparison, the inset shows the same data in linear scale. We use the same 



parameters as in Fig. 3(d) of the main paper: E 
ax = 0.05L and Or = 15°. 



lOhymL, hc/eB = L, Ri = L, R2 = 1.5L, Vo = 200hymL, Lt = 0.5L, 



